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ABSTRACT 

Sub milli-arcsecond astrometry and imaging of the black hole Sgr A* at the Galactic 
centre may become possible in the near future at infrared and sub-millimetre wave- 
lengths. Motivated by observations of short-term infrared and X-ray variability of Sgr 
A*, in a previous paper we computed the expected images and light curves, including 
polarization, associated with an compact emission region orbiting the central black 
hole. We extend this work, using a more realistic hot-spot model and including the ef- 
fects of opacity in the underlying accretion flow. We find that at infrared wavelengths 
the qualitative features identified by our earlier work are present, namely it is possible 
to extract the black hole mass and spin from spot images and light curves of the ob- 
served flux and polarization. At radio wavelengths, disk opacity produces significant 
departures from the infrared behaviour, but there are still generic signatures of the 
black hole properties. Detailed comparison of these results with future data can be 
used to test general relativity and to improve existing models for the accretion flow in 
the immediate vicinity of the black hole. 

Key words: black hole physics, Galaxy: centre, infrared: general, submillimetre, 
techniques: interferometric, polarization 



1 INTRODUCTION 

Testing general relativity in a regime where spacetime cur- 
vature is large remains one of the primary goals of obser- 
vational astronomy. Black holes provide a natural target 
in which to focus these efforts. Nonetheless, due to their 
extremely compact nature, an unambiguous signature of 
strong gravity has been elusive so far. 

A number of researchers have suggested that imaging 
an optically-thin background accretion flow could provide 



IBroderick fc Loebl J2005al) it was shown that for a typical 
accretion model, opacity will be significant at radio wave- 
lengths, substantially altering the shape and visibility of the 
black hole "shadow" cast by the black hole's photon capture 
cross-section. 

Multi-wavelength polarization observations of accret- 
ing black holes have also been proposed as a method of 
probing the black hole vicinity l|Broderick fc Loebl |2005a : 



a direct test of strong gravity (see, e.g 
iTakahash j|20oi 12004 : IBroderick fc Loebl 



Falcke et nUEjfjfi 



2005a). The angu- 



lar scale of the black hole in the Galactic centre (identified 
with the radio source Sgr A*) is 5-10 /ias, twice as large 
as the nuclear black hole in M87, and substantially larger 
than all other candidates. As a result, Sgr A* is the most 
promising candidate for high resolution imaging. Within the 
next decade it is expected that a Very Long Baseline Array 
(VLBA) of sub-millimetre telescopes will exist, providin; 
20 /ias resolution imaging capabilities iDoeleman fc Bowe: 
l2004lMivoshi et aljEoM) . Thus, theoretical efforts to com- 
pute the images of realistic accretion flows is warranted. In 
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IConnors et alJl980l : lLaor et al For optically-thin but 

geometrically-thick accretion flows, as expected in low lu- 
minosity galactic nuclei such as ou r Galactic centre, it was 
shown bv lBroderick fc Loebl i2005al) that comparisons of the 
polarization near and above the optically thick/thin transi- 
tion frequency is indicative of the black hole spin. 

However, observations of near-infrared (NIR) and X- 
ray flaring of Sgr A* have implied t hat the inner regions 
of the accretion flow are nonuniform llBaganoff et alJl200ll: 
[ Genzel et all 120031: IGoldwurm et al] 120031: lEckart etaU 
12004 IGhez et al]l2004) . The time scale of the variability, 
~ 10 min, is comparable to the period of the innermost 
stable circular orbit (ISCO), and thus suggestive of an or- 
biting hot spot. Therefore, it is likely that modeling of 
Sgr A* images wi ll requir e the inclusion of variability. In 
IBroderick fc Loebl {2005b) 

it was shown that the images 
and light curves associated with a hot spot could be used 
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to measure the mass and spin of the black hole. There are 
additional plans to use the Phase Referenced Imaging and 
Astrometry (PRIMA) instrument, currently under construc- 
tion, at the Very Large Telescope (V LT) to obtain sub-mill i- 
arcsecond astrometry of these flares dPaumard et al]|2005f) . 

Here we i m prove upon the calculations of 
iBroderick fc Loeb! 



(2005b) by making use of a more 
realistic hot spot model consisting of a localized population 
of non-thermal electrons (which are produced, for example, 
by a magnetic reconnection flare). In addition, we include 
the effects of disk opacity for a typical radiatively inefficient 
accretion flow model. The opacity systematically decreases 
the magnification and polarization fraction, smoothes out 
the polarization angle light curve, circularizes the motion 
of the centroid, and symmetrizes the time averaged spot 
images. 

Section [5] presents a summary of the computational 
methods used and the models for the accretion disk and hot 
spot. Sections |3| and contain the light curves and centroid 
paths for NIR and sub-millimetre frequencies, respectively. 
The effects of opacity are highlighted in section |3] Finally, 
section |S| summarizes our conclusions. 

In what follows, the metric signature is taken to be 
( — h ++), and geometrised unites are used (G — c — 1). 



a{M) 


n° (cm 3 ) 


T° (JO 


n 5U( cm 3 ) 


Pnth 





3 X 10 7 


1.7 x 10 11 


8 x 10 4 


-2.9 


0.5 


3 x 10 7 


1.4 x 10 11 


5 x 10 4 


-2.8 


0.998 


1 x 10 7 


1.5 x 10 11 


1 x 10 5 


-2.8 



Table 1. The accretion flow parameters associated with equa- 
tions H] 

terms of these, the radiative transfer equation is 



dA 



■- j u - a v N v 



(3) 



where j v and a v are emissivity and absorption coefficients 
(to see how these are related to their standard d efinitions, 
see IBroderick fc Blandfor3l2004l lLindauistlll96rt . The two 
synchrotron emission components consi dered here involve 
populations of electrons with a ther mal llYuan et alJl20o3) 
and power-law Jjones fc Q'Delll977D distributions. Because 
we are primarily concerned with emission at sub-millimetre 
or shorter wavelengths, Faraday rotation and conversion are 
unlikely to be important, and are thus neglected. 



2.2 Disk Models 



2 COMPUTATIONAL METHODS 

2.1 Ray Tracing & Radiative Transfer 

The method by which the light rays were produced and 
th e radiative transfer pe r formed is discussed at length 
in IBroderick fc Blandfordl i2003|) (for which tracing null 
geode sies is a limiting case) and IBroderick fc Blandfordl 
(2004) (sec also Broderick 200]J), respectively. As such, only 
a brief summary is presented here. 

Null geodesies are constructed by integrating the equa- 
tions 



da^ 
~dT 



/(r)F 



and 



dfc M 
dA 



-fir) 



1 dk v k v 

2 dx» 



(1) 

where the partial differentiation is taken holding the covari- 
ant components of the wave four- vector, constant. The 
function f(r) is arbitrary (corresponding to the freedom in- 
herent in the affine parameterization, A) and chosen to be 



f(r) = r* 



r 



(where th is the horizon radius) in order to regularize the 
affine parameter near the horizon. IBroderick fc Blandfordl 
( 2003) have explicitly demonstrated that this procedure does 
produce the null geodesies. 

Polarized radiative transfer in curved spacetime is most 
easily performed by integrating the Boltzmann equation 
jLindauistl Il965 IBroderick! l2005h . In this case, it is the 
photon distribution function N v oc I^/v 3 that is evolved. 
In the case of polarized radiative transfer, it is possible 
to define covariant analo gues of the Stokes paramet ers, 
N v = (N V ,N$,N?,N¥) IBroderick fc BlandfordlEooi) . In 



W e use the same thre e accret ion flow models as described 
in [ Broderick fc Loebl i2005bl) . Motivated by lYuan et alJ 
(2003), who showed that the vertically-averaged electron 
density and temperature are approximately power-laws in 
radius, we write the thermal electron density, n e , tempera- 
ture T e , and the non-thermal electron density n ntn as 



jj) exp(-2 /2p 



< M 



"nth = "-nth \ exp(-Z /2p ) . 



(4) 



where p is the cylindrical radius relative to the black hole 
spin axis, and p ntn is the radial power-law index for the non- 
thermal electrons. The assumed constants are listed in Table 
In all models, the non-thermal electrons have a spectral 
index of 1.25 and a minimum Lorentz factor of 100. The 
magnetic field strength is set to be a fixed fraction (30%) 
of equipartition relative to the protons. As suggested by re- 
cent general-relativistic magnet ohydrodynamic simula tions, 
the field is taken to be toroidal jDe Villiers et al]|2003T) . The 
accreting gas is assumed to be in free fall inside of the ISCO 
(6M, 4.233M and 1.237M for a = 0, 0.5 and 0.998, respec- 
tively), and in Keplerian rotation otherwise. In all cases the 
disk angular momentum was aligned with the spin of the 
black hole. 



2.3 Hot Spot Model 

The hot spot is modeled by an overdensity of non-thermal 
electrons centered at a point orbiting at the Keplerian ve- 
locity with a Gaussian profile as measured in the comoving 
frame: 



n s 



o 

n s exp 



v) 2 



2RI 



(5) 
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Figure 1. The spectral flux for the non-rotating (a = 0, solid), 
moderately rotating (a = 0.5, long-dashed), and maximally rotat- 
ing (a = 0.998, short-dashed) Sgr A* disk models as viewed from 
45° degrees above the equatorial plane. The thick lines correspond 
to the quiescent disk emission and the thin lines correspond to 
the peak flare flux for a spot located at an orbital ra dius of 6M. 
The d ata points are taken from the compilation by lYuan et alJ 
(2001). 




Figure 2. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the aver- 
age background subtracted total flux in the optically-thin regime 
as functions of time for spot circular orbits at 6M (solid), 8M 
(long-dash) and lOAf (short-dash) around a Schwarzschild black 
hole, viewed from 45° above the equatorial plane. The top panel 
shows the polarization angle, ifi. The time axis is set so that a 
single orbital period of the 10M case is shown. For a black hole 
mass of 4 X 10 6 M (as in Sgr A*), the time unit is M = 20s. 



where Ar M = r M — rg is the displacement from the spot 
centre (located at r-g), iig is the spot four- velocity, ng is 
the spot central number density, and Rs is a measure of the 
spot radius. So tha t our re sults may be easily compared with 
IBroderick fc Loebl J2005bl) . we used R s = 1.5M. The spot 
central number density was chosen to roughly reproduce the 
observed peak NIR flare flux. As suggested by coincident 
NIR-X-ray observa tions, the hot spot spectral index was 
assumed to be 1.3 fcckart et all 12004'). and the minimum 
Lorentz factor was assumed to be 100. Such a situation may 
arise, e.g., in the case of reconnection event similar to a so- 
lar flare. The spectrum of the disk and our canonical flare 
model (at an orbital radius of 6M, viewed from 45° above 
the equatorial plane) are shown in Figure Below 100 GHz 
the disk photosphere extends beyond the hot spot, quench- 
ing its contribution to the over all flux. Above 100 GHz the 
hot spot becomes increasingly visible, reaching a maximum 
near 500 GHz. 

For the Galactic centre, the orbital time scale near the 
ISCO are between 5 min and 30 min for orbits around max- 
imally and non-rotating black holes, respectively. The via- 
bility of a hot-spot model for the flaring in Sgr A* requires 
that the synchrotron cooling time scale be large in compari- 
son to the orbital period. Since the emission at a frequency v 
is dominated by electrons with Lorentz factor 7 e = \J v/vb, 
where vb is the classical electron cyclotron frequency, the 
cooling time of a flare as observed at v is 

2 f v \-i/2 / B \~ 3/2 



For the frequencies considered here this ranges from many 
hours (in the sub-millimeter regime) to tens of minutes (in 
the NIR). Hence, the highest frequencies at which a sub- 
stantial portion of the orbit is likely to be visible are those 
in the NIR. 

In order to concentrate upon the hot spot emission, all 
following plots are of background subtracted quantities, i.e., 
the quiescent disk emission has been removed. 



3 NEAR INFRARED 

As shown in IBroderick fc Loebl d2005al) . above ~ 10 3 GHz, 
the accretion flow is optically thin everywhere. As a result, 
at NIR frequencies the emission from the hot spot is unadul- 
terated. For this reason we discuss the light curves and cen- 
troid paths in the NIR first. Note that due to the fact that 
the spectra above ~ 10 3 GHz has constant spectral index, 
the results of this section apply to all wavelengths above 
the optically thick/thin transitio n. These may be dire ctly 
compared with the results of IBroderick fc Loebl J2005bl) . 

3.1 Light Curves 

As in IBroderick fc Loebl (l2005bl) . the primary distinction be- 
tween hot spots located at different radii is the time scale of 
variability in the total flux. This can be seen explicitly in the 
lower panel of Figure |21 in which the magnification of the hot 
spot emission is shown as a function of time. This is evident 
in the polarized flux variability as well (middle panel). It 
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Figure 3. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the aver- 
age background subtracted total flux in the optically-thin regime 
as functions of time for a spot orbit at 6M around a Schwarzschild 
black hole viewed from 0° (solid), 22.5° (dash-dot), 45° (long- 
dash) and 67.5° (short-dash) above the equatorial plane. The po- 
larization angle, ip, is shown in the top panel. For a black hole 
mass of 4 X 10 6 M (as in Sgr A*), the time unit is M = 20s. 



Figure 4. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the aver- 
age background subtracted total flux in the optically-thin regime 
as functions of time for a spot orbit at 6M around a Kerr black 
hole with a = (solid), a = 0.5 (long-dash) and a = 0.998 (short- 
dash) viewed from 45° above the equatorial plane. The polariza- 
tion angle, if), is shown in the top panel. The time axis is set so 
that a single orbital period of the a = case is shown. For a black 
hole mass of 4 X 10 6 M (as in Sgr A*), the time unit is M = 20 s. 



should be noted that the polarization variabi l ity will be dif- 
ferent from that found in iBroderick fc Loebl l)2005b|) owing 
to the difference in the polarized emission models (the one 
utilized here being the more physically motivated) . Nonethe- 
less, the generic features of a primary minimum caused by 
the development of the secondary hot-spot image (immedi- 
ately preceding maximum unpolarized magnification) and a 
weaker secondary minimum caused by the development of a 
tertiary image (following the maximum unpolarized magni- 
fication by 60M) exist. In this case the polarization angle is 
punctuated by rapid rotations (~ 90°) at the primary min- 
imum, and uneventful otherwise, rotating as expected for a 
toroidal field. 

Characteristically, the peak polarized flux follows the 
peak magnification. In the canonical case of an spot orbital 
radius of 6M viewed from 45° above the orbital plane, this 
corresponds to a delay of approximately 5 min in the context 
of the Galactic centre. 

Again, variations in viewing inclin ation, produces 
similar effects as those reported in IBroderick fe Loebl 
J2005bl) . the primary distinction being the total magnifi- 
cation. There are structural changes in the polarized flux, 
indicative of the promin ence of the primary and se condary 
hot-spot images. Unlike IBroderick fe Loebl (|2005b), in the 
edge-on view ($ = 0°), strong gravitational lensing produces 
a polarization component orthogonal to dominant polariza- 
tion, leading to a rapid rotation in polarization angle preced- 
ing peak magnification by approximately a quarter-phase. 

The manifestation of bl ack hole spin in the light c urves 
is shown in Figure H As in IBroderick fc LoebJ J2005bl) . the 
differences among the magnification light curves is small in 
comparison to those associated with variations of the orbital 



parameters. While the secondary minimum of the polariza- 
tion does vary significantly with spin, in the total flux this 
is difficult to observe. 

However, the primary distinction between different 
spins is likely to be the variations in the radius of the ISCO, 
and thus the typical periods. Figure^] shows the magnifica- 
tion and polarization light curves for a spot viewed from 45° 
above the equatorial plane located at the prograde ISCO 
for various dimensionless spin parameters (normalized by 
M), namely 6M, 4.233M and 1.237M for a = 0, 0.5 and 
0.998, respectively. Because the orbital timescales vary by 
nearly an order of magnitude between a — and a = 0.998, 
these are plotted as functions of orbital phase. Rapid rota- 
tion substantially reduces the maximum polarized flux, pri- 
marily due to the enhanced gravitational lensing associated 
with the compactness of th e orbit. However, as mentioned in 
IBroderick fc Loebl J2005bj) . and alluded to above, the most 
significant discriminator is likely to be the event time scale. 

3.2 Centroid Paths 

The PRIMA instrument at the VLT is expected to provide 
sub-milli-arcsecond astrometry, thus enabling high resolu- 
tion measurements of the the infrared image centroid for 
Sgr A* and in turn constrain the accretion flow and black 
hole parameters. Because the location of the image centroids 
will be dominated by the brightest features, subtracting the 
background accretion flow emission is necessary to isolate 
the location of the hot spot during a flare. However, intrin- 
sic variability in the source and/or the presence of multiple 
hot spots may introduce substantial systematic errors into 
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Figure 5. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the aver- 
age background subtracted total flux in the optically-thin regime 
as functions of phase for spot orbits at the prograde ISCO around 
a Kerr black hole with a = (solid), a = 0.5 (long-dash) and 
a = 0.998 (short-dash) viewed from 45° above the equatorial 
plane. The polarization angle, if), is shown in the top panel. The 
time axis is set so that a single orbital period of the a = case is 
shown. For a black hole mass of 4 X 10 6 M (as in Sgr A*), the 
time unit is M = 20 s. 



Figure 6. The paths of the background subtracted intensity cen- 
troid in the optically-thin regime for circular spot orbits around a 
Schwarzschild black hole viewed from 45° above the orbital plane 
with radii 6M (solid), 8M (long-dash) and 10M (short dash). For 
reference, a circle inclined at 45° is also shown by the dotted lines 
for each radius. Axes are labeled in units of M (corresponding to 
an angular scale of roughly 5 ^tas for Sgr A*). 



10 



this process. Nonetheless, in order to clearly remove the un- 
certainty associated with the accretion disk model, in com- 
puting all of the centroid positions below we have used the 
background subtracted images. These may be generated using 
the centroids of the total emission if the quiescent centroid 
position and intensity are known: 



— F(t)X(t) - FbqXbg 



Xs(t) 



F(t) - F BG 



(7) 



where Xs{t) is the location of the spot centroid, X(t) and 
F(t) are the location of the total image centroid and the 
total observed flux, and Xbg and Fbg are the background 
(disk) image centroid and flux. 

While gravitational lensing and relativistic aberration 
distort the centroid path, it is nevertheless possible to de- 
termine the orbital radius (Figure|SJ and inclination (Figure 
Combined with the orbital period at an orbital radius rs 
(in units of M), 



2ttM 



,3/2 



(8) 



(where positive a corresponds to prograde orbits) deter- 
mined, e.g., by the magnification light curve, this provides 
a direct measurement of the spin. 

Figures |H| and [5] present the effects of black hole spin 
upon the centroid paths. In Figure|H]the hot spot is located 
at 6Af and viewed from 45° above the orbital plane. In this 
case a moderate increase in the orbital velocity with spin 
leads to an enhancement of the relativistic beaming, shifting 
the centroid path towards the approaching side (left). In 
Figure [5] the hot spot is located at the prograde ISCO for 
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Figure 7. The paths of the background subtracted intensity cen- 
troid in the optically-thin regime for circular spot orbits around a 
Schwarzschild black hole with radius 6Af viewed from 0° (solid), 
22.5° (dash-dot), 45° (long-dash) and 67.5° (short-dash) above 
the orbital plane. For reference, a circles inclined at each angle 
are shown by the dotted lines. Axes are labeled in units of M 
(corresponding to an angular scale of roughly 5 ^tas for Sgr A*). 
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Figure 8. The paths of the background subtracted intensity cen- 
troid in the optically-thin regime for circular spot orbits with ra- 
dius 6M in the equatorial plane around a Kerr black hole viewed 
from 45° above the orbital plane for a = (solid), a = 0.5 (long- 
dash) and a = 0.998 (short-dash). For reference, a circle inclined 
at 45° is also shown by the dotted line. Axes are labeled in units 
of M (corresponding to an angular scale of roughly 5 fias for Sgr 
A*). 



10 




10 

-10 



-5 



10 



Figure 9. The paths of the background subtracted intensity cen- 
troid in the optically-thin regime for circular spot orbits around 
a Kerr black hole viewed from 45° above the orbital plane at 
the prograde ISCO for a = (solid), a = 0.5 (long-dash) and 
a = 0.998 (short-dash). For reference, a circle inclined at 45° is 
also shown by the dotted lines for spin. Axes are labeled in units 
of M (corresponding to an angular scale of roughly 5 fias for Sgr 
A*). 



each black hole spin, again viewed from 45° above the orbital 
plane. Here the centroid path at high spin, where the ISCO 
is close to the horizon, is substantially offset. However, in the 
case of moderate spin, the generic dependence upon radius 
seen in Figure HJ is dominant. 



4 GENERIC EFFECTS OF OPACITY 

While in the NIR the disk opacity to synchrotron self- 
absorption is negligible, this is not the case at sub-millimeter 
wavelengths, at which the accret ion flow is only beginning t o 
be optically thin. As discussed in lBroderick fc Loebl (l2005al) . 
the accretion flow opacity is not symmetric between the two 
sides of the disk, and is enhanced by the Doppler effect on 
the approaching side of the disk. As a result, the light curves 
and centroid paths presented in the previous section for the 
NIR can be substantially modified for frequencies near the 
optically thick/thin transition. For the purpose of highlight- 
ing the generic effects of opacity, and the utility of high- 
resolution multi-wavelength flare observations, we present 
next a comparison of the light curves and centroid paths for 
NIR and radio frequencies. 



4.1 Time Averaged Images 

The effects of asymmetric opacity are readily apparent in the 
orbit-averaged images, shown in Figure ITT71 for the canonical 
hot-spot model for the three disk models discussed in Sec- 
tion Ell as observed at 230 GHz, 350 GHz and in the NIR. 
In the NIR both the direct and secondary images are clearly 
visible, corresponding to the larger and smaller rings, respec- 
tively. Less visible are the tertiary images, which produce the 
very thin ring directly inside that associated with the sec- 
ondary images. In addition, as expected from the Doppler 
shift and relativistic beaming, the approaching portion of 
the hot-spot orbit (left side) appears considerably brighter 
than the receding portion of the spot orbit (right side). At 
radio frequencies the secondary image is substantially less 
visible and the brightness contrast between the approaching 
and receding portions of the hot-spot orbit decreases with 
decreasing frequency. This is a direct result of the increased 
opacity in the approaching portio n of the disk ( also on the 
left side). As mentioned in iBroderick fc Loebl i2005bl) . in 
the limit of a thermal spectrum the orbit-averaged images 
are symmetric, and thus it is not surprising that increased 
opacity results in more symmetric images. 

Also shown in Figure 1101 is the degree and orientation 
of the polarized flux. In all cases this is dominated by the 
approaching portion of the orbit. In the NIR, a second po- 
larization component is present due to the secondary im- 
age. This is substantially suppressed at radio wavelengths, 
and thus the strong swings in polarization angle are not ex- 
pected in the radio. Finally, the degree of polarization is 
also expected to decrease as a result of the decrease in the 
brightness asymmetry mentioned before. 



4.2 Light Curves 

Figures I11H13I shown the magnification and polarized flux 
light curves in the radio and NIR for a hot-spot orbiting an 
a = 0, 0.5 and 0.998 black hole, respectively. The primary 
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Figure 10. Comparison of the orbit-averaged, disk-subtracted images of a spot for two radio frequencies (at which opacity is important) 
and in the infrared (at which the disk is transparent). For reference, the photon-capture cross section for a Schwarzschild black hole is 
shown by the dashed white line. 



distinction between the different observing wavelengths is in 
the degree of polarization, decreasing by roughly an order 
of magnitude from the NIR to 230 GHz. In addition, as ex- 
pected from the images, the polarization angle light curve 
is progressively smoothed as the observing wavelength in- 
creases. It should also be noted that the secondary minimum 
in the polarization, also resulting from the development of 
tertiary images, is significantly less visible in the radio. 



4.3 Centroid Paths 



The centroid paths for different observing frequencies are 
compared in Figures [TUllbl As the observing frequency de- 
creases the centroid path becomes larger and more symmet- 
ric. As a result, the orbital parameters are more easily deter- 
mined at radio wavelengths. This implies that simultaneous 
high-resolution radio and NIR observations of the centroid 
are capable of identifying black hole spin. 
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Figure 11. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time for a 
spot orbit at 6M around a Schwarzschild black hole viewed from 
45° above the equatorial plane as observed in the infrared (solid), 
350 GHz (long-dash) and 230 GHz (short-dash). The polarization 
angle, tp, is shown in the top panel. The time axis is set so that a 
single orbital period of the a = case is shown. For a black hole 
mass of 4 X 10 6 Mq (as in Sgr A*), the time unit is M = 20s. 
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Figure 12. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time for a 
spot orbit at 6M around a Kerr black hole (a = 0.5) viewed from 
45° above the equatorial plane as observed in the infrared (solid), 
350 GHz (long-dash) and 230 GHz (short-dash). The polarization 
angle, tp, is shown in the top panel. The time axis is set so that a 
single orbital period of the a = case is shown. For a black hole 
mass of 4 X 10 6 M (as in Sgr A*), the time unit is M = 20s. 
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Figure 13. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time for a 
spot orbit at 6M around a Kerr black hole (a = 0.998) viewed 
from 45° above the equatorial plane as observed in the infrared 
(solid), 350 GHz (long-dash) and 230 GHz (short-dash). The po- 
larization angle, is shown in the top panel. The time axis is 
set so that a single orbital period of the a = case is shown. For 
a black hole mass of 4 X 10 6 Mq (as in Sgr A*), the time unit is 
M = 20 s. 
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Figure 14. The paths of the background subtracted intensity 
centroid for circular spot orbits around a Schwarzschild black hole 
viewed from 45° above the orbital plane at 6M as observed in the 
infrared (solid), 350 GHz (long-dash) and 230 GHz (short-dash). 
For reference, a circle inclined at 45° is also shown by the dotted 
lines for each radius. Axes are labeled in units of M (correspond- 
ing to an angular scale of roughly 5/ias for Sgr A*). 
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Figure 15. The paths of the background subtracted intensity 
centroid for circular spot orbits around a Kerr black hole (a = 
0.5) viewed from 45° above the orbital plane at 6M as observed 
in the infrared (solid), 350 GHz (long-dash) and 230 GHz (short- 
dash). For reference, a circle inclined at 45° is also shown by 
the dotted lines for each radius. Axes are labeled in units of M 
(corresponding to an angular scale of roughly 5 ^tas for Sgr A*). 



10 



-10 



350 GHz 




10 



10 



Figure 16. The paths of the background subtracted intensity 
centroid for circular spot orbits around a Kerr black hole (a = 
0.998) viewed from 45° above the orbital plane at QM as observed 
in the infrared (solid), 350 GHz (long-dash) and 230 GHz (short- 
dash). For reference, a circle inclined at 45° is also shown by 
the dotted lines for each radius. Axes are labeled in units of M 
(corresponding to an angular scale of roughly 5 ^tas for Sgr A*). 



Figure 17. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time at 
350 GHz for spot orbits at 6M (solid), 8M (long-dash) and 10 M 
(short-dash) around a Schwarzschild black hole, viewed from 45° 
above the equatorial plane. The polarization angle, ip, is shown in 
the top panel. The time axis is set so that a single orbital period 
of the 10M case is shown. For a black hole mass of 4 X 10 6 Mq 
(as in Sgr A*), the time unit is M = 20s. 



5 RADIO 

Anticipating high-resolution sub-millimetre observations of 
the Galactic centre, we present the magnification and po- 
larized flux light curves and image centroids at radio fre- 
quencies. As discussed in the previous section, these may be 
expected to be quantitatively different from those presented 
in Section [3] 



5.1 Light Curves 

Figure HTI shows that at 350 GHz the absorption has little ef- 
fect upon the magnification light curves for moderate orbital 
inclinations (compare to Figure|5J. However, the peak polar- 
ized flux is reduced by more than 70% and observed earlier 
than in the NIR. While the polarization angle continues to 
exhibit a rapid rotation associated with the development of 
a secondary image, it is significantly smoothed. 

This is especially apparent in Figure 1181 which shows 
the light curves for a number of viewing inclinations, 
(compare to Figure In the case where ■& = 0°, the emer- 
gence of an orthogonal polarization is suppressed, leading to 
a constant polarization angle. For small #, the peak magnifi- 
cation is also marginally reduced, as may be expected from 
the larger column density for such viewing angles. Never- 
theless, the flare time-scale and peak magnification are still 
indicative of the orbital radius and inclination. Furthermore, 
the presence of a flare in the polarized flux is generic. 

Figures 1191 and 1201 show the magnification and polar- 
ized flux light curves for various black hole spins as seen 
at 350 GHz and 230 GHz, respectively (compare to Figure 
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Figure 18. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time at 
350 GHz for a spot orbit at 6M around a Schwarzschild black 
hole viewed from 0° (solid), 22.5° (dash-dot), 45° (long-dash) 
and 67.5° (short-dash) above the equatorial plane. The polariza- 
tion angle, ip, is shown in the top panel. For a black hole mass of 
4 X 10 6 M (as in Sgr A*), the time unit is M = 20 s. 



Figure 20. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time at 
230 GHz for a spot orbit at 6M around a Kerr black hole with 
a = (solid), a = 0.5 (long-dash) and a = 0.998 (short-dash) 
viewed from 45° above the equatorial plane. The polarization an- 
gle, ip, is shown in the top panel. The time axis is set so that a 
single orbital period of the a = case is shown. For a black hole 
mass of 4 X 10 6 Mq (as in Sgr A*), the time unit is M = 20s. 
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Figure 19. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of time at 
350 GHz for a spot orbit at 6M around a Kerr black hole with 
a = (solid), a = 0.5 (long-dash) and a = 0.998 (short-dash) 
viewed from 45° above the equatorial plane. The polarization an- 
gle, ip, is shown in the top panel. The time axis is set so that a 
single orbital period of the a = case is shown. For a black hole 
mass of 4 X 10 6 Mq (as in Sgr A*), the time unit is M = 20s. 



0J. Typically, there is little difference among the magnifica- 
tion light curves. However, in both cases the polarized flux 
clearly discriminates between high and moderate/low black 
hole spins, varying by roughly a factor of two in both cases. 
This is likely due to the lower thermal electron density per- 
missible in the high spin model employed (see Table 0. 

As in the optically-thin regime described in Section r3.ll 
the primary difference between high and low spin black holes 
is likely to be the flaring timescales. Similarly to Figure |S| 
Figure VTW shows the magnification and polarized flux light 
curves as functions of orbital phase. As in the NIR, the po- 
larization is more sensitive to the black hole spin than the 
degree of magnification. The opacity associated with the 
hot-spot orbiting the maximally rotating black hole is signif- 
icantly larger due primarily to the rapid disk velocity at this 
position. As a result, this model was significantly more ab- 
sorbed than the others. Combined with strong gravitational 
lensing, this effect produces a variation in the flux that is 
unpolarized and small in comparison with the quiescent disk 
emission. 

5.2 Centroid Paths 

Next, we consider the paths of the background subtracted im- 
age centroids. While in the optically thin regime described 
in Section lo\2l the disk and hot-spot emission are completely 
disentangled, in this case the disk opacity can substantially 
modify the image centroids. In the frequency regime of the 
optically thick/thin transition, opacity is most significant on 
the approaching side of the accretion flow, near where the 
hot-spot is the brightest. This has the effect of reducing the 
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Figure 21. The background subtracted total and polarized flux 
(bottom and middle panels, respectively) normalized by the av- 
erage background subtracted total flux as functions of phase at 
350 GHz for spot orbits at the prograde ISCO around a Kerr 
black hole with a = (solid), a = 0.5 (long-dash) and a = 0.998 
(short-dash) viewed from 45° above the equatorial plane. The po- 
larization angle, if), is shown in the top panel. The time axis is 
set so that a single orbital period of the a = case is shown. For 
a black hole mass of 4 X 10 6 Mq (as in Sgr A*), the time unit is 
M = 20 s. 



Figure 22. The paths of the background subtracted inten- 
sity centroid at 350 GHz for circular spot orbits around a 
Schwarzschild black hole viewed from 45° above the orbital plane 
with radii 6M (solid), 8M (long-dash) and 10M (short dash). For 
reference, a circle inclined at 45° is also shown by the dotted lines 
for each radius. Axes are labeled in units of M (corresponding to 
an angular scale of roughly 5 fia,s for Sgr A*). 



10 



brightness contrast between the approaching and receding 
portions of the hot-spot orbit. Since the brightest regions 
dominate the location of the image centroid, this has signif- 
icant implications for the centroid paths. 

Figure l2~2l shows the centroid paths for hot-spots orbit- 
ing at three radii, viewed from 45° above the orbit al plane. 
This is similar to Figure 12 in iBroderick fc Loebl <2005bl) . 
primarily due to the reduced brightness contrast mentioned 
above. The orbits are noticeably larger than those in Figure 
|S| though the major axis distance is similar. 

Various viewing inclinations are shown in Figure l2*o1 (cf.. 
Figure |7J . As in Figure 1221 the centroid paths are larger 
than those found in the NIR. For low #, this substantially 
modifies the shape of the centroid path. Nonetheless, the 
nearest region of the orbit (bottom) is close to the unlensed 
orbit positions (shown by the dotted lines), and thus the 
semi-minor axis is still diagnostic of the orbital inclination. 

Figures l?4l and 1231 show the centroid paths for various 
black hole spins as viewed at 350 GHz and 230 GHz, respec- 
tively (cf., Figure [SJ. In all cases the paths are larger than 
in the NIR and, unlike the NIR, are larger than expected 
in the absence of gravitational and opacity effects (shown 
by the dotted line). In all cases, the semi-minor axis again 
appears to be a good measure of the orbital inclination. 

The centroid paths for orbits at the ISCO of the three 
black hole spins we considered are shown in Figure 1261 As 
in the NIR, the inner orbits are substantially modified by 
strong lensing, and in this case, opacity. However, as men- 
tioned previously, the flux variation associated with the 
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Figure 23. The paths of the background subtracted inten- 
sity centroid at 350 GHz for circular spot orbits around a 
Schwarzschild black hole with radius 6Af viewed from 0° (solid), 
22.5° (dash-dot), 45° (long-dash) and 67.5° (short-dash) above 
the orbital plane. For reference, a circles inclined at each angle 
are shown by the dotted lines. Axes are labeled in units of M 
(corresponding to an angular scale of roughly 5 ^tas for Sgr A*). 
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Figure 24. The paths of the background subtracted intensity 
centroid at 350 GHz for circular spot orbits with radius 6M in 
the equatorial plane around a Kerr black hole viewed from 45° 
above the orbital plane for a = (solid), a = 0.5 (long-dash) and 
a = 0.998 (short-dash). For reference, a circle inclined at 45° is 
also shown by the dotted line. Axes are labeled in units of M 
(corresponding to an angular scale of roughly 5 (J,as for Sgr A*). 
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Figure 25. The paths of the background subtracted intensity 
centroid at 230 GHz for circular spot orbits with radius 6M in 
the equatorial plane around a Kerr black hole viewed from 45° 
above the orbital plane for a = (solid), a = 0.5 (long-dash) and 
a = 0.998 (short-dash). For reference, a circle inclined at 45° is 
also shown by the dotted line. Axes are labeled in units of M 
(corresponding to an angular scale of roughly 5 ^tas for Sgr A*). 
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Figure 26. The paths of the background subtracted intensity 
centroid at 350 GHz for circular spot orbits around a Kerr black 
hole viewed from 45° above the orbital plane at the progradc 
ISCO for a = (solid), a = 0.5 (long-dash) and a = 0.998 (short- 
dash). For reference, a circle inclined at 45° is also shown by the 
dotted lines for spin. Axes are labeled in units of M (correspond- 
ing to an angular scale of roughly 5 /ias for Sgr A*). 

high-spin case is small, and thus it is not clear that the 
centroid path will be measurable in this case. 



6 CONCLUSIONS 

The photon trapping radius of the black hole in the center of 
the Milky- Way galaxy has the largest apparent sizes among 
all known black holes, occupying an anglular scale of tens 
of micro-arcseconds on the sky. This scale will be within 
reach of forthcoming observatories, such as a sub-millimetre 
VLBA (for imaging) , or the PRIMA instrument on the VLT 
(for monitoring shifts in the infrared image centroid during 
flares). Imaging hot spots in the immediate vicinity of an ac- 
creting black hole provides a new method for testing general 
relativity and measuring the black hole mass and spin. 

The image and centroid shifts of a compact population 
of non-thermal electrons may be used to constrain the mass 
and spin of the central black hole. At radio wavelengths such 
a spot would persist for many orbitally periods, ultimately 
being sheared into a ring, while in the NIR such a spot would 
cool rapidly, surviving approximately a single orbit. 

At the sub- millimetre frequencies currently proposed for 
the VLBA imaging of the Galactic centre, the opacity of 
the underlying steady accretion flow has a significant effect 
upon the images and centroid motions. Nonetheless, it re- 
mains possible to extract the black hole mass and spin from 
these. Furthermore, combined with the NIR data (at which 
the accretion flow opacity is negligble) this will provide a 
method by which the accretion flow may be isolated. 

While polarization is significant in both the radio and 
NIR regimes, the enhanced opacity at radio wavelengths 
reduces the polarization signatures of strong lensing that 
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are most noticeable for the polarization angle and polarized 
flux in the N1R. Immediately prior to maximum magnifi- 
cation, rapid changes in the orientation of the polarization 
take place in the N1R. These are systematically smoothed in 
the radio, though the polarization does rotate as expected 
for a toroidal magnetic field. 

We have used simple models for the underlying ac- 
cretion flow and the hot spot. Future work can improve 
upon our results using three-dimensional, fully relativistic, 
magneto-hydrodynamic simulations of the accretion flow, in 
which hot spots (or flares) arise naturally as magnetic re- 
connection events and/or at shocks. 
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